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Abstract. In multiple time-scale (singularly perturbed) dynamical systems, 
canards are counterintuitive solutions that evolve along both attracting and 
repelling invariant manifolds. In two dimensions, canards result in periodic 
oscillations whose amplitude and period grow in a highly nonlinear way: they are 
slowly varying with respect to a control parameter, except for an exponentially 
small range of values where they grow extremely rapidly. This sudden growth, 
called a canard explosion, has been encountered in many applications ranging 
from chemistry to neuronal dynamics, aerospace engineering and ecology. Canards 
were initially studied using nonstandard analysis, and later the same results were 
proved by standard techniques such as matched asymptotics, invariant manifold 
theory and parameter blow-up. More recently, canard-like behaviour has been 
linked to surfaces of discontinuity in piecewise-smooth dynamical systems. 

This paper provides a new perspective on the canard phenomenon by showing 
that the nonstandard analysis of canard explosions can be recast into the 
framework of piecewise-smooth dynamical systems. An exponential coordinate 
scaling is applied to a singularly perturbed system of ordinary differential 
equations. The scaling acts as a lens that resolves dynamics across all time- 
scales. The changes of local curvature that are responsible for canard explosions 
are then analyzed. Regions where different time-scales dominate are separated 
by hypersurfaccs, and these are pinched together to obtain a piecewise-smooth 
system, in which curvature changes manifest as discontinuity-induced bifurcations. 
The method is used to classify canards in arbitrary dimensions, and to derive the 
parameter values over which canards form either small cycles (canards without 
head) or large cycles (canards with head). 
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1. Introduction 

This paper revisits the geometry of dynamical systems exhibiting relaxation 
oscillations. Extreme changes of flow curvature brought about by the presence of 
different time-scales are studied. After discussing previous approaches, we introduce 
new methods to characterize these changes. We show that a correspondence exists 
between certain singularities of current interest in two different kinds of dynamical 
system: those that are singularly perturbed j and those that are piecewise-smooth 

A system with two time-scales can typically contain regions of fast dynamics, 
evolving rapidly towards lower dimensional manifolds where regimes of slow dynamics 
take over. The stability of a slow manifold depends on it being hyberbolically 
attractive in its normal direction (9) [TUl [2Ql [21] ■ The occurrence of non- hyperbolic 
points leads to the separation of the fast and slow dynamics becoming indistinct. 
This is well-known to generate such nontrivial behaviours as relaxation oscillations 
and canard explosions [30 . The term "canard" alludes to the duck-like shape of 
cycles that wrap around critical manifolds of a typically cubic form (look ahead to 
figure [2] panel (2)), but also to the deception they play on the mathematician by their 
sudden appearance and disappearance. The explosion is a cascade of cycle growth 
that occurs in a parameter range that is exponentially small in the ratio between 
time-scales (a small parameter e), implying that it is very difficult to observe either in 
practical or numerical situations, hence it could easily be construed as a discontinuous 
event. How small the time-scale ratio (e) must be for this to occur is addressed in a 
related paper [8]. In this paper we show how such a discontinuous description can be 
derived, and still preserve sufficient geometry to capture the canard phenomenon in 
detail. 

An extreme view of a slow-fast dynamical system can be taken by expressing it 
in terms of ordinary differential equations containing discontinuities. A vector field 
defines evolution throughout phase space, but its value can jump discontinuously 
across a surface called a switching manifold. In the systems we will consider, orbits can 
slide along the switching manifold in a manner that approximates the slow dynamics, 
while fast dynamics takes place outside the switching manifold. The conventional 
description of dynamics in piecewise-smooth vector fields was laid down by Filippov 
We will rederive Filippov's convention by different means, defining it as a 
piecewise-smooth union of vector fields of the fast and slow subsystems, combined 
with a map that 'pinches' phase space in the neighbourhood of the critical manifold. 
The pinching creates a hypersurface, the switching manifold, and the slow dynamics 
slides around inside it. 

The pinched system compresses the cascade of cycles in a canard explosion into a 
single instant. We make a precise connection between this pinching formalism and 
the methods of nonstandard analysis applied to singular perturbation theory [1]. 
We begin by combining pinching with the nonstandard concept of an exponential 
microscope, and we show that approximations obtained by both methods are 
equivalent. Furthermore, an observation from [1 allows us to make successive 
approximations that separate the canard explosion into two cascades - one with small 
cycles ( "canard without head" ) and one with large cycles ( "canard with head" ) - each 
of which is instantaneous in the pinched system. 

We demonstrate the pinching method on a classic example of a slow-fast system, 
namely the two-dimensional van der Pol oscillator [28]. We then apply pinching to a 
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system with two fast variables and one slow, which allows for oscillations on the fast 
time-scale observed as "bursting" oscillations. This is the classical Hindmarsh-Rosc 
burster [15] . where a spike-adding mechanism takes place and is organised by saddle- 
type canard cycles. The same mechanism can be found in other three-dimensional 
bursters (such as the Morris-Lecar equations [27]) and has been recently studied from 
the standpoint of invariant manifolds, in particular, slow manifolds of saddle-type and 
their (un)stable manifolds [13] . We demonstrate that, by using pinching microscope, 
one can understand the transition from canard cycles to spiking periodic attractors, 
that is, limit cycles with one spike. 

The systems we analyse can all be expressed in the form 

ex + xf'(x) =g(x,y,z), (1) 

where e is a small positive constant, and /' is the x-derivative of a smooth function 
/, which is a cubic polynomial in x. The function g depends on x, and on a variable 

y = ex + f(x), (2) 

and may depend on a third dynamic variable z. We can then write ([T]) as a system of 
ordinary differential equations, 

sx = y- f(x), (3a) 
y =g{x,y,z), (3b) 

from which it is clear that, for e <C 1, the dynamics will separate into a slow time-scale 
t and a fast time-scale t/e. In cases where we have a third equation i = k(x, y, z), the 
dynamics of z may have either a fast or slow time-scale; the example we consider in 
section l4~3l falls into the first category. 

The paper is arranged as follows. In section [2] we introduce some classic canard 
theory based on singularly perturbed systems, including the nonstandard analysis 
concept of an exponential microscope. Then in section [3] a description of canards 
in terms of pinched (or piecewise-smooth) vector fields is introduced, and applied to 
examples in section [4] In sections [4.11 and [4721 we consider supercritical and subcritical 
cases of the van der Pol oscillator, the latter representing an abstract form of the 
Fitzugh-Nagumo equations, and in section 14.31 we apply the method to describe the 
birth of spiking via canard cycles in a three dimensional Hindmarsh-Rose model. In 
section[S]we apply pinching to derive a piecewise-smooth description of non-hyperbolic 
singularities and some of their related canards. Some concluding remarks are made in 
section El 



2. The classic canard example 

2.1. The van der Pol equations in the Lienard plane 

Results in this section derive from the seminal paper on canards published by the 
French group from Strasbourg in the early 1980's [1|. Thoughout, they consider the 
equations of the van der Pol oscillator, which can be written as 

ex = y — ]-x 3 + x, (4a) 

y =q-x, (4b) 

in the Lienard plane (x, y) [23], with constant forcing q and singular perturbation 
parameter < e -C 1. 
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The system ^ exhibits a supercritical Hopf bifurcation as q passes through 
q = ±1, that is, a stable focus- type equilibrium loses stability by creation of a stable 
periodic orbit. The location of the small parameter e in Q distinguishes the nullcline 
x = 0, which we label So and call the critical manifold. It consists of attracting 
and repelling branches, S a and S r , adjoined at the turning points (or folds) of So, at 
[x, y) = ±(1, -f ), such that S = S a US r . 

Let us first consider how the critical manifold organises the dynamics of the 
system for e = 0, illustrated in figure [T] Setting e = in ^ gives a differential- 
algebraic system consisting of the differential equation y = q — x on the slow time- 
scale, constrained to So by the algebraic equation y — x 3 /3 — x. Hence So is the phase 
space of this limiting problem, called the "slow subsystem" (or sometimes "reduced 
system" ) . 




Figure 1. Fast and slow subsystems of the van der Pol equation. Double arrows 
indicate a fast motion. The fold points of the critical manifold So correspond to 
the bifurcation points of the fast dynamics; they separate the attracting sheets 
S a of the critical manifold from the repelling sheet S r . 

To study the fast subsystem one introduces a fast time-scale r = t/e. Denoting 
differentiation with respect to r by a prime, Q becomes 

x = y — —x 3 + x, (5a) 

y' = £{q - x). (5b) 

When e = 0, the slow variable y has its dynamics frozen, that is, y remains constant 
and can be considered as a parameter. This yields a family of one- variable differential 
equations on the fast time-scale, parametrised by y and usually referred to as the "fast 
subsystem" (or sometimes "layer problem"). The full system with e = is piecewise- 
smooth, with a switch occurring between the fast and slow subsystems. From (|5ap 
we observe that the (y-parameterized) equilibria of the fast subsystem belong to none 
other than the critical manifold So, being stable on S a and unstable on S r (indicated 
by the filled and unfilled disks in figure Q}. 

The aim in this paper, as in singular perturbation theory, is to build from this 
insight on the limiting case e = to study what happens when e is small but nonzero. 
Here we recount the main results, summarized in figure [21 which shows a simulation 
of (0J| for e = 5 X 10~ 3 . If an initial condition is taken a distance of 0(1) away 
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from the critical manifold Sq, the solution of the van der Pol equation first evolves 
on a fast time-scale, almost horizontally, until it reaches an e-neighbourhood of one 
of the attracting sheets S a . This first epoch is well captured by the fast subsystem 
of horizontal fibres in figure [1] Once in the neighbourhood of S a , where x sa 0, the 
slow subsystem begins to dominate the dynamics and the orbit moves slowly, close to 
S a . The sole event that can end this slow epoch is a loss of normal stability of the 
attracting sheet S a of the critical manifold. This corresponds to a bifurcation point 
of the fast subsystem, generically a saddle-node bifurcation, and takes place at either 
of the folds in figure [TJ 



Canard with head Maximal canard Canard without head 




0.95 0.97 0.99 q 1.01 

Figure 2. Canard explosion in the van der Pol system. The central panel (a) 
shows the bifurcation diagram of this system as a function of the parameter q 
(I/2-norm || ■ || of the limit cycle on the upper branch, and equilibria along the 
lower branch, unstable on the dashed part to the left of HB). The quasi- vertical 
part of the branch of limit cycles corresponds to the canard explosion. Five dots 
(other than the Hopf point HB) are marked along the branch: © is a stable 
equilibrium solution and (2) to © are periodic solutions, of which (2) to © are 
canards and © is a relaxation oscillation. The five solutions are represented in 
the state space (x,y) in the outer panels, together with the cubic nullcline So- 

A simple stability analysis of ([4]) reveals a unique fixed point where x = y = 0, 
which lies on So for all q. For \q\ > 1 the fixed point is stable and lies on S a , for 
\q\ < 1 it is unstable and lies on S r . When q = ±1, the system undergoes a Hopf 
bifurcation which takes place at the fold of So, creating a branch of stable periodic 
orbits that displays the so-called "canard explosion" . A surprising property of the 
canards becomes clear when one plots orbits for various values of q together with 
the critical manifold So in the Lienard plane. As shown in figure [H a periodic orbit 
follows the attracting sheet S a of the critical manifold down to the fold point, then, 
counterintuitively, instead of being ejected at the fold along a fast fibre, the orbit 
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sticks to So along its repelling sheet S r for an 0(1) time. 

The family of canard cycles are distinguished as being "with head" or "without 
head" as in figure [2J The distinction depends on whether an orbit leaves the e- 
neighbourhood of S r along a fast fibre moving to the left (with head) or to the right 
(without head) as depicted. The intervening case occurs at a parameter value q = qo 
and is called the maximal canard, shown in panel 3 of figure [5J The maximal canard 
follows the unstable region for the greatest time possible by traveling the whole length 
of S r between the folds. The value of go can be estimated numerically or calculated 
analytically using asympotic expansions in e (e.g. [5]), and in section f4. II we estimate 
it geometrically. 

To understand why the canards stick to the repelling critical manifold S r , it is 
usual to consider them from the viewpoint of invariant manifolds. In the singular 
limit e = 0, Sq is an invariant manifold. The question of the persistence of 
invariant manifolds for e > was addressed in the 1970's in the work of Hirsh, Pugh 
and Shub |16] and, in the context of slow- fast dynamical systems, in the work of 
Fenichel (9] [10] ; see also [20] [21] . Fenichel proved that compact normally hyperbolic 
subsets of an invariant manifold, such as So in the present case, persist as locally 
invariant manifolds S e (generally non-unique) for every small enough e > 0, and that 
they are smooth e-perturbations of the unperturbed manifold. Normal hyperbolicity 
can be lost, however, as happens generically at (isolated) fold points in planar systems, 
where attracting and repelling sheets of invariant manifolds intersect. Therefore, in 
the case of the one-dimensional cubic critical manifold So, one can apply Fenichel 
theory everywhere except the fold points. 

The key to the canard explosion is the fact that a small parameter change can alter 
the arrangement of the attracting and repelling branches S® and S^. of the invariant 
manifold S e . Specifically, the size of a canard cycle is determined by whether the 
arrangement of S*° and S £ forces it to curve inside the critical manifold So (figure |3](a)) 
to form a canard without head, or outside it (figure EHc)) to form a canard with head. 
Since the two invariant manifolds exchange their position between figures [3] (a) and 

(a) (b) (c) 



q = 0.995 t' = 0.05 = 0.9934909325 £ = 0.05 o = 0.99 £ = 0.05 




Figure 3. Invariant manifolds in the neighbourhood of the fold during the canard 
explosion. Showing the attracting and repelling invariant manifolds 5° and S^, 
and the critical manifold So = S a US r , for different values of q. In (a), 5° curves 
inside allowing for a small canard cycle; in (b), and SJ coincide; in (c), the 
two manifolds have exchanged their positions and the resulting periodic attractor 
is a large canard cycle. 

(b) as q varies, there must exist a parameter value q = qo at which they pass through 
each other. Then, as the manifolds are themselves orbits, they must coincide as shown 
in figure |3fb). The canard that passes from 5" to S^ will remain on the repelling sheet 
S r for the longest distance possible, creating the maximal canard shown in panel @ 
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of figure [5J The strong (exponential) repulsion in the normal direction from S£ causes 
the transition between small (without head) and large (with head) canard cycles to 
occur within an interval of q that is of order e~ c l e for some c > 0. This implies that 
it is not trivial to resolve different canard cycles, because their repelling segments 
are all contained in an exponentially small thickness around S r . In order to obtain 
quantitative as well as qualitative information about the canards along the explosive 
bifurcation branch, changes of coordinate were introduced in [1] that magnify the 
geometry close to S r , such that different canards within the same subfamily (i.e. with 
or without head) can be separated. The crucial point is that this magnification must be 
exponentially strong in e in order to distinguish between orbits that are exponentially 
close to each other, as we show in the following section. 



2.2. Exponential microscope 

To understand relaxation oscillations in more detail we must separate the canard 
cycles, figureEl from the critical manifold, Sq. We follow [Tj by defining a new variable, 
v = x, which has the effect of stretching the phase portrait by a factor of 1/e in the 
direction transverse to Sq. In the (x,v) plane, (d]) becomes 

x = v, (6a) 

ev = q — x — (x 2 — l)v. (6b) 

The dynamics in these coordinates is somewhat changed in appearance, however it 
is clear that there is a one-to-one correspondence between orbits of the van der Pol 
system in the Lienard plane (x,y) and in the phase plane (x,v). In these coordinates 
a; is a slow variable compared to v, so the slow dynamics is organised around a critical 
manifold v = 0, (in place of So which is now the v = axis), given by 

q — x 

S v = {(x 1 v) : v = j(x)}, where j(x) = -= -. (7) 

x z — 1 

The change to coordinates (x, v) implements a magnification transverse to So by 
a factor e. Although this separates the canard cycles from So (as we shall see in 
section f4. II) . it is not strong enough to separate canard cycles from each other, and in 
particular cannot resolve the maximal canard. To this end, Benoit et al. introduced 
an exponential microscope. For this stronger magnification we must do better than 
zooming in on Sq. As observed in [T], the maximal canard can be approximated near 
the fold by S v when q = 1, which by (O is given the hyberbola v — 70(2;) = — jj—. 
Thus the required 'microscope' variable is defined as 

w = (v - l0 (x)) l£ \ (8) 

employing the notation 

a:M = \x\ p sgn(x) (9) 

which will recur throughout the paper. As was found in [1] (and we rederive in 
section l4~Tj) using it? = ^ • and letting 5 = {q — qo)^- E \ we now have 

x = w [1/£] +70, (10a) 



X - qp 
7o 



(10b) 



recalling that qo is defined as the parameter value for which the canard is maximal 
(figure [2] panel @). 
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The term w (S /iv) is a perturbation of the 6 = (maximal canard) system. 
It is small for \w\ 3> S, where the perturbation is negligible, but in the strip \w\ <C S 
this term dominates and the vector field is almost vertical. The illustration of this in 
figure 2] is based on the original analysis in [T], which can be summarized as follows. 




Figure 4. The van der Pol system under the exponential microscope: (a) q < qo, 
{8 < 0), (b) q = q , (8 = 0), (c) q > qo, (8 > 0). 



In the parameter regime of the canard explosion, asymptotic analysis reveals that 
5 = q — qo = exp(— c/e) for a constant c > 0, guaranteeing that 8 <C 1. At the 
parameter value q = qo, when the canard cycle is maximal, (jTUJ) reduces to 

x=w [1/£] , (11a) 

w = w . (lib) 

7o 

Noticing that qo lies in an e-neighbourhood of unity, and that canard orbits are 
exponentially close to the hyperbola v = Jo(x) near the fold, we can approximate 
the system as 

x w - — j— , (12a) 
1 + x 

-w(x-l)(x + l). (12b) 

System (TT12)) is integrable, and provides the following approximation for the solutions 
of system (TlT]) . 

— + Y~T~ X ) , (13) 

where K is a constant. In figure [4] we show the phase portrait given by (|13p outside 
the strip \w\ > 5. Inside \w\ <C S, the phase portrait is approximated by vertical fibres 
with a direction given by the sign of S. When we discuss the method of pinching 
in section 14.11 we will show how keen an insight this was, with computed plots that 
reveal just how sharply the vector field changes between the strip's exterior, \w\ > S, 
and its interior, \w\ < 6. 



3. Pinched dynamical systems and sliding 

In this section we introduce a piecewise-smooth formulation of multiple time-scale 
systems, obtained by a method called pinching. We will use this later to complement 
and extend the analysis in section [2. 2 1 The method of pinching has its origins in [5], 
where it was used to overcome the lack of uniqueness of orbits in Filippov's description 
[11] of piecewise-smooth dynamical systems. 

Consider a smooth vector field F : 1" i— >• R™ , and the dynamical system 

x = F(x). (14) 
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Let there exist a smooth scalar function h : U H > M on an open region W C R n , such 
that V/i ^ everywhere on Let /i have two level sets labeled 

E± = {x 6 K" : h(x) = ±a} , (15) 

for some a > 0, bounding a strip \h\ < cr which we call the pinch zone. (For motivation, 
consider the pinch zone to be the strip of quasi- vertical flow, \w\ < S, in figure@|. We 
shall denote by h the Lie derivative of h along the flow of (TH)) , 

h = x- — h(x) = F Vh. (16) 
ax 

We wish now to approximate the solutions of x = F in the pinch zone by some rule 
that: (i) associates each point on E + with a point on and (ii) replaces the vector 
field between them by a suitable approximation. 

In the local region U, we can define a surjective map from the pinch zone 
|/i(x)| < cr, to a switching manifold £ C K n_1 . We do this by choosing fibres x = pj (A) 
that connect the two boundaries £+ and £_, given by 

{x = p c (A) : p e (+l) e £+, P€ (-l) e E_, -1< A < 1} , (17) 

so that each fibre is associated with (or projected onto) a unique point (eE. Two 
such fibres are illustrated in figure [5f a). The vector field on £ is then the set 

7- e = {F(p e (A)):-KA<l}. (18) 

For simplicity we assume that, for any £, the set J-j contains at most one vector that 
is tangent to a level set of h; we call this a sliding vector, and denote it as 

F 5 (0 - {F(p c (A)) : M P5 (A)) = 0, |A| < l} . (19) 

We label the hypersurface h = as S. For F$ to be unique one must have that, for 
a given £, the chord Pf(A) is never tangent to S, hence we choose p such that the 
derivative dh/dX is non- vanishing, which we can write as the condition 

|^( Pc (A)) = • V(F(x) • VMx)) ^ 0. (20) 

Moreover, if J-(£) contains a sliding vector, F^^), then we describe £ as a sliding 
point, otherwise it is a crossing point. Over any open region of sliding points, the 
sliding vectors form an n — 1 dimensional vector field, and thus define a dynamical 
system on £ given by 

e=Ps(0- (2i) 

Orbits of the pinched dynamical system are defined as the concatenation of: 

(i) solutions of equation (TH)) in the regions h(x) > a and h(x) < —a, and 

(ii) solutions of equation (|2"Tj) on £. 

The following theorem follows directly from (TT9"|) - ([2"T]) . 

Theorem 3.1. (a) The boundary between crossing and sliding lies where F is 
tangent to £ + or £_. (b) Zeros ofF inside the pinch zone correspond to zeros of the 
sliding vector field Fs(£), which we call sliding equilibria. 

Taking h as a local coordinate, a vector field for the pinched system is obtained 
by making a piecewise-smooth coordinate transformation h i— >■ h, where 

h = h- crsgn(/i) for \h\ > a, h = for \h\ < a. (22) 

We say that the sliding region is: 
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Figure 5. Pinching, and the origin of sliding, (a) A smooth vector field F, and 
a pinch zone bounded by S±, which are joined by fibres p^(A). The set-valued 
vector field is shown on fibres corresponding to a sliding point £1 and a 

crossing point £2- Sliding is determined by the positioning of a tangency surface 
S (where h = 0), on which we take the sliding vector F$. Above and below 
the pinch zone the vector field is approximated by F±. (b) The pinched system 
around a switching manifold S, obtained by pinching together S_|_ and S— . A 
tangency of F+ to £ bounds sliding (right of S) and crossing (left of 5) on £. 

• stable if /i(pj(+l)) > > h(p(— 1;£)) (orbits flow locally towards S), and 

• unstable if h(p(— 1;£)) > > /i(pf(+l)) (orbits flow locally away from E). 

The intersection of S with £± is the locus of points where F is tangent to E+ or 
£_. It is this intersection that partitions E into regions of sliding or crossing, so the 
domain of sliding is given by S (1 {x : \h(x)\ < a} . 

Figure [5] illustrates how the dynamics of a smooth system (a) is approximated, 
after pinching, by a piecewise-smooth system (b). Orbits that cross both boundaries 
£± of the pinch zone in the same direction in (a), are approximated by orbits crossing 
the switching manifold E in (b), while orbits that cross the boundaries in opposite 
directions in (a) give rise to sliding on the switching manifold in (b). Also illustrated 
is the content of theorem 3.1, showing that the boundary between crossing and sliding 
is a tangency between the flow and the pinch zone boundary in (a), corresponding 
to tangency between the flow and the switching manifold in (b), and showing that 
an equilibrium inside the pinch zone in (a) is approximated by an equilibrium in the 
sliding vector field in (b). 

3.1. Filippov's sliding vector field 

Equation (f2Tj) requires knowledge of the exact form of the vector field F inside the 
pinch zone, but in many circumstances a further approximation is useful. Recalling 
that J-j is the set of values of F(x) along the fibre p^(A), we can approximate it with 
a convex set JF^ t which is an interpolation between the values of F at the endpoints 
Pf(±l). Expressing these endpoint values as F±(£) = F(p^(±l)), we thus have 

F € « ^ ee | i±^F + (0 + i^F_ (0 : -1< A < 1 J . (23) 

This is precisely the convex combination of two vector fields F± used by Filippov 
[TTj to define sliding on a switching manifold. A sliding vector, which lies on the 
hypersurface S, is a vector F e T that satisfies F • V/i = 0. Assuming that 
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V/i(pj(+l)) ~ V/i(pf(— 1)), this occurs at A = X s , where 

_ (F_+F + )-Vfr 
5 (F_-F+)-Vft ' 

and therefore the sliding vector field is given by 

- (F_ ■ V/QF+ - (F + ■ Vfr)F_ 

Fs = (F_-F + ).Vft ■ (24) 

This is formally identical to Fillipov's sliding vector field, and provides us with a 
piecewise-smooth dynamical system 

r F+, if fc(x) > a, 
±=\ F S , if |fc(x)| <*, (25) 
[ F , if h(x) < -a. 

How well (1231) approximates the smooth system (|14p clearly depends on how the pinch 
zone | h \ < a is chosen. For a slow-fast system such as the van der Pol oscillator 
(|4al) with small e, we can choose the pinch zone to be a neighbourhood of the critical 
manifold So- Importantly, the interesting quantities we calculate in later sections 
depend on the geometry of S, which itself depends on the choice of 'pinch function' 
h, but not on the width of the pinch zone a. 

We will now apply pinching to take a fresh look at some familiar examples of the 
canard phenomenon, before using it to derive a classification of canards in section [5l 

4. Applications of pinching in 2 and 3 dimensions 

4-1- Pinching and sliding in the van der Pol oscillator 

In this section we use the pinching formalism to analyse canards in the van der Pol 
oscillator. In figure [5](a) we illustrate a relaxation oscillation simulated from (0| for 
e = 0.04. The critical manifold Sq is given by h — 0, where 

h(x, y) = y - ^x 3 + x. (26) 



(a) (b) (c) (d) 




Figure 6. Simulation of {4{ for e = 0.04 and q = 0.9, showing the same relaxation 
oscillation (bold) in four coordinate systems: (a) the Lienard plane (x,y), (b) 
the plane (x,h), (c) the exponential microscope (x, V), (d) the piecewise-smooth 
system (x, V). The critical manifold, consisting of attracting and repelling 
branches S a and S r , is straightened to give (b), a pinch zone \V\ < e e emerges 
around it in (c), which is pinched to a switching manifold V = in (d). 

If we vary the parameter q from the situation in figure 151(a), there will be no observable 
change in the segment of periodic orbit lying on the attracting righthand branch of the 
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critical manifold S a , prior to its sudden disappearance in a canard explosion (recall 
figure [5]). Thus we attempt to resolve different orbits by proceeding along similar lines 
to the exponential magnification in section [521 but we choose to magnify the Lienard 
variable y directly instead of using the phase variable v = x; this has the benefit that 
x remains the fast variable. First, we simply straighten out the critical manifold So 
by taking ft as a variable, giving the system, 

ex = h, (27a) 

h = q- x- (x 2 - 1)-. (27b) 

e 

The relaxation oscillation from figure [^a) is shown in the (x, h) plane in figure HJb). 
The critical manifold Sq is now the horizontal axis h = 0. 

We now attempt to resolve canard orbits by magnifying the neighbourhood of So . 
To do this we introduce a variable 

V = ee sgn(h)h e : (28) 

giving the system 

ex =V [1/£ \ (29a) 

V = V ((1 - x 2 ) + V [ - 1/£] (q - x)ej , (29b) 

in which the same relaxation oscillation is plotted in figure [5tc). The critical 
manifold Sq remains on the horizontal axis, V = 0, and is evidently no longer 
a local attractor/repeller, but is surrounded by a region, approximately \V\ < e £ 
(corresponding to \h\ < e), where the flow is quasi- vertical. To consider the dynamics 
relative to this strip, we find the nullcline V — for V ^ 0, given from (|29bl) by 



S v = {(x,V)eR 2 : V={e{x-q)/{l-x 2 )Y i }. (30) 
This is shown in figure [7] for different values of q. The curve S v becomes singular 
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Figure 7. Exponential microscope. Phase portrait of the van der Pol oscillator 
in the (x, V) magnification, for: (a) q < 1, (b) q = qv = 1, (c) q > 1. The dashed 
curve is the nullcline V = 0. 

when (dV/dx,dV/dV) = (0,0), forming the cross seen in the lower right region of 
figure [Tfb), at a parameter q — q v . From (a) to (c) the nullcline Sv thus undergoes 
a bifurcation [6] in the form of a saddle transition. Solving (dV/dx, dV/dV) = 
gives the coordinates of the singular point as (x, V) — (q v , — (e/2) e ), and substituting 
this into ([29b]) gives q v = 1, which therefore coincides with the Hopf bifurcation in 
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the original system. The relaxation oscillation in figure [TJa) disappears in a canard 
explosion between (a) and (b). 

We now derive a piecewise-smooth approximation from (|29[) by following section[31 
to arrive at figure [HJd). It is natural to choose a pinch zone \V\ < e £ (equivalently 
\h\ < e), inside which the vector field is near vertical. We will pinch along the In- 
direction (setting p in section [3] as Px(A) = (x,\e £ )). The vertical fibration of the 
pinch zone by the flow implies that this is a good approximation to the dynamics. This 
is equivalent to introducing a new coordinate as in (|22|) . given by V — V — e £ sgn(V). 
From section [21 sliding occurs at x values for which the nullcline S v (V = 0) lies in 




-2 



Figure 8. Relaxation oscillation in: (a) the (x, V) magnification, (b) the 
piecewise-smooth system obtained by pinching the region \V\ < e E in (a). Plotted 
for the same parameter values as figure Oa). So is the V = axis, the dashed 
curve is the nullcline V = corresponding to h = 0. 

the pinch zone. From (|30[) . the sliding region is therefore given by 

e £ <e £ => \x-q\ < |1 - x 2 \, (31) 

on the switching manifold V = 0. The sliding vector field has an equilibrium at 
x = q. The boundaries of sliding, where the vector field is tangent to the switching 
manifold, are therefore where |a? — g| = |1 — x 2 \. The tangencies to the upper and 
lower side of the manifold are labeled T\ and respectively in figure [9] The critical 
manifold So has become the switching manifold, with stable sliding corresponding 
to the attracting manifold S a , and unstable sliding corresponding to the repelling 
manifold S r , as shown. 

It is not necessary to calculate the sliding vector field explicitly (though it is 
easily found from as we can infer it qualitatively from two facts: firstly that the 

stability of the sliding equilibrium (stable for \q\ < 1 and unstable for \q\ > 1), and 
secondly that, from (f24| . the sliding vector field is equal to the upper vector field at 
Ti and to the lower vector field at T2. 

Inspection of the piecewise-smooth vector field in figure |U(a) makes it clear that 
a relaxation oscillation occurs because an orbit, sliding towards T2 from the right, will 
detach from the manifold at T2, pass through the manifold near x = — 1, and then be 
re-injected into the sliding region from above. The bifurcation of the curve S v at q = 1 
causes the tangency points to collide in figure |9jb). Then the crossing region (dotted) 
between T\ and T2 shrinks to nothing, and a canard orbit exists that passes from 
the stable to unstable sliding regions (from S a to S r ). Every point on the unstable 
sliding region S r belongs to a family of orbits that travel along S r for some time, 



W\ = 



x — q 



l-x 2 
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Figure 9. Sketch of the bifurcation responsible for the canard explosion: (a) 
relaxation oscillation at q < 1, (b) the simple canards at q = 1, (c) stable sliding 
equilibrium at q > 1. The upper and lower vector fields are tangent to the 
switching manifold at Ti and T2 respectively. These bound the sliding regions 
(horizontal lines) and crossing regions (dotted horizontal lines). The black spot 
is the sliding equilibrium, unstable in (a) and stable in (c). S a and S r label 
the stable and unstable sliding regions, approximating stable and unstable slow 
manifolds of system (I29D . 



before departing into either the upper or lower half planes and being reinjected into 
S a ; every one of these is a periodic orbit and constitutes an instantaneous 'snapshot' 
of the entire canard explosion. When the tangency points separate in figure HJc) , 
orbits in the righthand sliding region are attracted to the sliding equilibrium and no 
periodic orbits are possible. Thus the bifurcation that takes place in figure [9] describes 
the sudden disappearance of the relaxation oscillation, via the instantaneous existence 
of a family of canard cycles at q = 1. In section [5] we will classify this as the simple 
canard case of a catastrophic sliding bifurcation. 

Let us now show that the same method can split the instant q = 1 into two families 
of canards: those in figure [DJb) that detach from the upper side of the switching 
manifold (without head) and those that detach from the lower side (with head). 

Whereas in figures [6] (a) and (b) the cycle (bold) appears to lie on an attracting 
branch S a of the critical manifold as it approaches the fold, we see from figure U^sl), 
in which S a lies on the axis V = 0, that the cycle lies instead on the curve V = 0. 
Let us therefore approximate the canard by the curve S v when 5 = 1. Solving h = 
(equivalent to V — 0) for q — 1, from (|27b[) , confirms the result that led to (O, that 
near the knee of So at (x,y) — (1, —2/3), the canard is given by 

h « ejo(x) = -£——■ (32) 
1 + x 

To magnify the phase portrait around this curve in the /i-direction we use a different 
exponential scaling, introducing the variable 

W = (h-e l0 (x)) [e] , (33) 
and obtaining the system 

x = -W [ - ] +70 (a), (34a) 
s 

W= \W\ x -i ([w 1 * 1 +670(3:)] [(1 - x 2 ) - £70(2;)] + (g - a:)e).(34b) 

The relaxation oscillation in the (x, W) plane is illustrated in figure [TOta). The focus, 
which in the (x, V) plane resided inside the pinch zone, has now been lifted outside. 
Compare this to figure Hta), which is restricted to — 1 < x < 2 and a vertical range 
that excludes the focus; the variable W differs from the analysis of Benoit et al leading 
to ((5|) by a scaling W = e £ w. Notice that the vector field, while vertical in the region 
given approximately by \ W\ < e 2s , changes direction as it passes through the nullcline 
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W — 0. As q changes from (a) to (c) in figure [TU1 the nullcline segment inside the pinch 
zone sweeps to the right, changing the vertical vector field direction from downward 
to upward. 
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Phase portrait of the van der 



Second exponential microscope. 
Pol oscillator in the (x, W) magnification, for: (a) q < q w , (b) q = q w , (c) 
1 > q > q > W. The dashed curve is the nullcline V = 0. Compare the region 
x > to figure 3] 



The nullcline = is given from (|34bp by 
5 W = J (a, W) e M 2 



W = 



x — q 



1 - x 2 - £7o0) 



70 (z) 



(35) 



Similar to the bifurcation of the nullcline S v at g = 1, the curve becomes singular 
at some parameter q = q w , forming a cross as shown in figure HUT b) . This takes 
place when (dW/dx, dW/dW) = 0, the solution of which gives the coordinates of the 
singularity as 



(x,W) 



qw, 



3e 2 



2(1 + q w ) (q w (l + q w f - e) 




(36) 



and substituting these into W = gives the bifurcation parameter q — q w as the 
solution of 

e = (1 - q w )(l + q w f . (37) 

Similarly to the (x, V) system, we now define a pinch zone given by \W\ < e 2s in 
which the vector field is almost vertical. We pinch vertically, in this case by defining 
a new variable W = W — e 2s sgnW, and in the (x, W) plane we obtain a piecewise- 
smooth system as shown in figure QTJ We omit the details for calculating the sliding 
dynamics on the switching manifold W — 0, as the steps are similar to those for the 
(x, V) system, and because the necessary information can be inferred directly from 

The vector field is tangent to the switching manifold at the x-values where the 
curve (1551 enters the pinch region. We label them T\ and T 2 for the tangency to the 
upper and lower side of the switching manifold respectively. As illustrated in figure [T2l 
for q < q w , T\ lies to the left of T2, with stable sliding to the right of T2. This allows 
a periodic orbit to exist as simulated in figure ITTT b') and sketched in figure [TST a). As 
the nullcline W = bifurcates, at q = q w given by (|37[) . the two tangencies collide. A 
canard orbit then passes from the stable to unstable sliding regions, and since every 
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(a) (b) 
W=0 




Figure 11. Relaxation oscillation in: (a) the (x,W) magnification, (b) the 
piecewise-smooth system obtained by pinching the region \W\ < e 2e in (a). 
Plotted for the same parameter values as figure Oa) . The dotted curve is the 
nullcline W = 0. 



point in the unstable sliding region belongs not only to an orbit in the sliding region, 
but also to orbits departing into both the upper and lower half planes, every one of 
these is also a periodic orbit and constitutes an instantaneous 'snapshot' of a canard 
explosion. This explosion is not the one observed in the {x, V) plane, figure[9j because 
it terminates not in a stable focus, but in the small periodic orbit shown in figureJT^Ic). 
This, instead, is a cascade between the canard with head in figure [T27 a) and the canard 
without head in figure 112( c) . In section [3] we will classify this as the visible canard 
case of a catastrophic sliding bifurcation. 




Figure 12. Sketch of the bifurcation responsible for the transition between 
canards with and without head: (a) a large cycle corresponds to the canard 
with head at q = qw — 10~ 4 , (b) the visible canard at q = q w , (c) a small 
cycle corresponds to the canard without head at q = qw+ 10 — 4 . The upper and 
lower vector fields are tangent to the switching manifold at fold points T± and 
Ti respectively; these exchange exchange ordering in the bifurcation. S a and S r 
label the stable and unstable sliding regions, approximating stable and unstable 
slow manifolds of system J34D . 

Previously we assigned the value q = qo to the maximal canard, which marks the 
transition between canards with and without head. The bifurcation of the W = 
nullcline at q = q w signifies a change in curvature of the flow, from pushing orbits 
upwards to pushing them downwards relative to the maximal canard at W = 0. 
Therefore q w gives us a geometrical approximation for the value q = qo at which the 
maximal canard occurs, qo — q w , and hence from (j37[) we obtain a quartic expression 
for q : 

e = (1 - g )(l + go) 3 - (38) 
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Thus the bifurcation in figure [T2] controls the disappearance of a large cycle 
(figure H2T a)) to be replaced by a small cycle (figure [T2T c)). via an instantaneous 
canard cascade at q = q w . 

For q ss go we can use an observation of [I] that led to (|10b[) , namely that we can 
rewrite (|34a[) using the following result which is proved in the Appendix: 

Lemma 5.1 [lj The ordinary differential equation can be rewritten as 



1 w x-q q-q , , 

w w = — + wm £ - ^ 

Using this lemma, we can approximate for q — qo sufficiently small, 

x « 7q (x), (40a) 

W^W^— (40b) 
7o (x) 

Replacing the lcfthand side of p0b|) with W = x^- = 70 {x)W^^-, we can solve 
to find 

W(x) « W(0) exp (-q x - ^^x 2 + ^x 3 + (41) 

as first found in Q] and recounted in (fT3)) , except that here qo is included with the 
value (|38|) ) instead of approximation by unity. This gives an approximate expression 
for the trajectories of canards in the region \ W\ > e 2e . 

This can be improved using the pinched piecewise-smooth system, because the 
vector fields in the half planes W > and W < are independent. In particular we 
can make independent approximations in either half plane. We linearize about the 
equilibrium in the upper half plane, and expand for small W in the lower half plane 
(where there are no equilibria), and also expand about x = 1, to obtain 

x w - - + H(W)W [1/E] /e, (42a) 

W « 2W ((q -x) + H(W)(q - q ) (l - -W^ e ^j \ , (42b) 

where H(W) equals 1 if W > and if W < 0. This adds a correction to (|40| outside 
\W\ > e 2s that captures the focus in the upper half plane, and captures the exchange 
of tangency points responsible for the canard. Using the pinching method, we have 
thus derived a piecewise-smooth system, (|42|). that captures the geometry behind the 
canard explosion in figure [12] and extends the nonstandard analysis summarized in 
section [221 

Let us briefly review what information about the singularly perturbed system has 
been captured by pinching. The exchange of tangency points T\ and T2 in figures [§] 
and [12] enacts a discontinuity-induced bifurcation - bifurcations of the periodic orbits 
caused by interaction with the discontinuity at V = or W = 0. These cases are 
catastrophic forms of sliding bifurcations [18] . whose classification we will present 
in section [5] These sliding bifurcations in the (pinched) piecewise-smooth system 
correspond to bifurcations in the topology of the nullclines V = and W = in the 
smooth system, coinciding with the Hopf bifurcation and maximal canard respectively. 
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4-2. Pinching and sliding in the Fitzugh-Nagumo equations 

We now extend the analysis of the previous section, where an oscillator underwent a 
supercritical Hopf bifurcation, to a case involving a subcritical Hopf bifurcation. As a 
specific example we take an abstract form of the FitzHugh-Nagumo equations [12] used 
to model electrical excitations in nerve membranes. Written in terms of a membrane 
voltage potential v and recovery variable w, setting f(v) — /3v 3 + (1 + a)v 2 — av and 
g = g(v,w) affine, these take the form 

ev = /3v 3 + (1 + a)v 2 -av-w + I, (43a) 
w = bv — cw, (43b) 

where a,/?,/, are parameters and b, c, are constants. This is equivalent to ^ up to a 
simple coordinate translation, 

l + a I-w,(l + a\ .... 

x = v -^ y= ^- f {-w)> (44) 

giving Q with f(x) = ^x 3 — rx and g(x, y) — q — px — y, to obtain 

ex = y — -x 3 + rx, (45a) 
y =q-px-y. (45b) 

The critical manifold Sq of this system is the curve h — 0, where 

1 „ 

h(x,y) = y- -x +rx. (46) 

The form of So depends on r, and we fix r > so that Sq has two turning points 
(where f'(x) — 0), at (x,y) — ± A /r(l,— |). Throughout this section we will keep r 
and p fixed and vary q as a bifurcation parameter. 

The dynamics depends upon the number of equilibria of (|45[) . These occur at 
(x eq , y oq ) where x cq is a solution of ^x 3 q + (p — r)x cq — q — 0. The number of solutions 
depends on the discriminant A = (3q/2) 2 + (p — r) 3 . If A < there are 3 equilibria, 
and a stable cycle that can be destroyed by homoclinic connection to a saddle. In this 
section we restrict attention to the simpler case with one focus, A > 0, so we fix 

A = (3q/2) 2 + (p- r) 3 > 0. (47) 

The flow is illustrated in figure [13] A relaxation oscillation exists in (a). It is 
subsequently destroyed, but by a different process to that in section 14.11 There is 
a focus equilibrium which in (a) lies on the unstable branch S r of Sq. The focus 
changes stability in a subcritical Hopf bifurcation, creating an unstable cycle as shown 
(dotted) in figure fT5Tb) . The unstable cycle undergoes a canard explosion, and (b) 
captures an instant in the sequence of canards. In this system, the unstable canard 
is without head and grows in amplitude, while the stable canard is with head and 
shrinks in amplitude. In (c), the two cycles coincide at a parameter q = qa to be 
determined, and annihilate in a saddle- node bifurcation of cycles. At q = go, the two 
cycles coincide and we observe numerically that these are close to a maximal canard 
(dotted). Afterwards, in (d), the cycles are destroyed, leaving a stable focus on S a . 

We will not go into as much detail as the previous section, and instead seek only 
to derive the mechanism by which the cycles are destroyed in figure [T3T c), following the 
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Figure 13. Simulation of l|45|l for e = 0.04. Plots are with p = r = 1 and q 
values: (a) q - 0.02, (b) q - 6.045 X 10 -5 , (c) q , (d) q + 0.02. The cubic 
curve is So given by x = h = 0, with attracting branches S a and repelling branch 
S r . Stable (bold) and unstable (dotted) periodic orbits are shown around a focus 
(spot). The maximal canard parameter <yo is to be determined. 



same steps as in section I47TI We first straighten out So by transforming to coordinates 
(x, h), giving the system 

ex = h, (48a) 

h = r] p (x) + ^r]' s (x)h, (48b) 
in terms of the function 

i]x(x) = h(x, q - Xx) = ~^-x 3 + (r - X)x + q, (49) 

and its derivative rj' x {x) — r — A — x 2 . The nullcline h = is the curve h — e r y(x), 
where 

7 (a:) = -rtp(x)/rf e (x). (50) 

This nullcline bifurcates when dh/dx = dh/dh = 0. Solving these we find that the 
bifurcation takes place at coordinates 

z{e - p) 



(x, h) = {x h , h h ) = \ ±Vr - £, ^ J , (51) 

and substituting these coordinates back into h = we find the bifurcation parameter 
value to be 

q = qh = r^l + (p - r ) x h- (52) 

This coincides with the value of q at which an subcritical Hopf bifurcation takes place. 

We can magnify the strip \h\ < e by rescaling to a coordinate V = as in (|28p. 
obtaining 

ex = V [1/e \ (53a) 

V = \V\ (v'M + eW- 1/£l 77 p (x)) , (53b) 

then derive a pinched system to approximate it. The piecewise-smooth system 
obtained at this level is identical to figure HI that is, the unstable cycle belongs to 
the infinite number of coexisting periodic orbits at q = q^, none of which can be 
distinguished as either stable or unstable. As before, to distinguish between different 
canards we must magnify around the maximal canard. 
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The maximal canard, similarly to section 14.11 is found to lie on the curve given 
by the nullcline h = when q = q^. We therefore let 7o(x) be the function 7(2;) from 
(|50[) with q — qh, and approximate the maximal canard by h ~ £70(2;), where 
_ /s _ 'JyM - g = 1 / | 3p-2r-£ \ 

T]' e (x) 3 V x±y/r - ej 

The nullcline V = bifurcates at the same parameter q = qh as the nullcline h = 0. 

To magnify around this maximal canard curve, we rescale to a coordinate 
W = (h — £7o(x))[ £ l ) obtaining 

x = -W [1/e] + 70(a), (55a) 

£ 

- IVI^ 1 -* 1 ([^(x) - £ 7 ^(x)][W^ 1 + £ 7o (x)] + E7fc(a;)) , (55b) 

where 7o(a;) = [2x70(3;) + (p — r + x 2 )]/(r — £ — x 2 ). 

The nullcline W — bifurcates when dW/dx = dW/dW = 0, which we solve 
to find the bifurcation coordinates (x,W) = (iEo,Wo). For the same reasoning as in 
section 14.11 this bifurcation takes place approximately when the maximal canard is 
formed. We find that xq is a solution of 

^(x )-£/o>o)=0, (56) 
which can be solved numerically, or expanded as a Taylor series to find 

x ^ Xh + ^ { -s){&r~-7e) Xh + ° t {X0 ~ Xh) ^ (57) 

and, exactly, 

Wo = (-mixo) - x „ ( V \ iX0) ll( , ) ^ ■ (58) 

Substituting these into W = we find the bifurcation parameter q — q at which the 
maximal canard exists to be given by 
1 3 

Qa = 3^0 + ( r ~P)xo- (59) 

Simulations in this case are similar to figure [TTJ a), except that they reveal an unstable 
cycle around a stable focus as shown in figure rUT a) for q < <7q. At go in (b) , the 
tangencies T± and T2 of the upper and lower vector fields to the switching manifold 
collide, allowing a canard orbit to pass from the stable to unstable sliding regions, 
giving rise to a family of canard cycles. As for the supercritical van der Pol system, we 
have canards with head that depart from the underside of the unstable sliding region, 
and canards without head that depart from the upper side, but now these correspond 
to the cascades of stable and unstable canards, one of which is the maximal canard. 
Afterwards, for q > qo in (c), the piecewise-smooth system contains no mechanism by 
which a periodic orbit can exist, as all orbits are attracted towards the stable focus. 

Similarly to the supercritical van der Pol system, the bifurcation in the pinched 
system will be classified in section [5] as a visible canard case of catastrophic sliding 
bifurcation. We could go on to consider similar systems, such as replacing (|45bl) with 
y = q — px + y, in which a supercritical Hopf bifurcation is followed by a canard 
explosion that can be interrupted by a homoclinic connection to a saddle. In such 
a case, it might be appropriate to perform an exponential magnification around the 
stable manifold of the saddle to determine when it becomes part of a canard. Instead 
of considering further such extensions, let us now turn our attention to canards in 
three dimensions. 
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Figure 14. Sketch of the bifurcation responsible for the stable and unstable 
canard explosions in the subcritical van der Pol system: (a) a stable cycle (bold) 
with head and unstable cycle (dotted) without head at q = qa — 10 -4 , (b) a visible 
canard exists at the maximal canard parameter q = qo, (c) no periodic orbits exist 
for q = qo + 10 — 4 . The upper and lower vector fields are tangent to the switching 
manifold at fold points T\ and T2 respectively; these exchange exchange ordering 
in the bifurcation. 



4-3. Two fast variables in three dimensions: the Hindmarsh-Rose burster 

We now give an example of pinching in a three dimensional system. Unsurprisingly 
the analysis is more involved in three dimensions, and it is beyond the scope of this 
paper to investigate such a system in detail. We give only a brief introduction to 
the analysis that can be carried out, taking the Hindmarsh-Rose equations j!5| as an 
example. 

A slow-fast system in K 3 with two fast variables has a fast subsystem that is two- 
dimensional. This allows for the possiblity of periodic motion in the fast dynamics, 
which can induce bursting oscillations in the full system [17]. The Hindmarsh-Rose 
equations are a classic example of a simple bursting system, and can be written as 

ex = z — f(x) + 1 — y, (60a) 

y =s(x — xi) — y. (60b) 

ez = g{x) — z, (60c) 

where f(x) = ax 3 — bx 2 and g{x) = c— dx 2 . Here x represents a voltage variable and y 
a slowly varying current, giving a slow- fast system of the form ([3]), to which we add a 
fast variable z and applied current /, see [14] for more details. We take typical values 
a = 1, b = 3, c = 1, d = 5, and s = —1.618, that are known to be physically relevant. 

This model exhibits spiking behaviour. One can observe by numerical simulation 
that periodic bursting attractors are formed by an alternation between rapid phases, 
containing a number of spikes, and slow resting phases. Parameter-dependent families 
of periodic orbits are known to exist and display a spike-adding phenomenon. In the 
following we consider branches of periodic attractors obtained upon variation of /, but 
a similar spike- adding phenomenon exists when I is fixed and e is varied. At low I- 
values, the unique stable equilibrium of the system undergoes a Hopf bifurcation where 
a branch of stable periodic attractors is born at the point labelled y sp in figure [T5] 
The Hopf point is located very close to one fold of the w-nullcline and the subsequent 
periodic orbit follows the unstable (middle) branch of this cubic nullcline. Hence these 
periodic orbits are of canard type and the phenomenon is reminiscent of a canard 
explosion. For increasing J-values, the associated periodic orbit stays close to the 
unstable branch for a longer time, until eventually it remains close long enough to 
reach the upper fold point. This is consistent with the notion of a maximal canard. 
Shortly after this, as / increases, the associated periodic orbit starts to develop a fast 
oscillation - a spike. Hence the spike's appearance is suggested to be canard-induced, 
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and it is found that the process repeats itself again and again, as / increases, each 
time adding an extra fast spike, as shown in hgure [T51 




Figure 15. Periodic attractors of the Hindmarsh-Rose model 11601 1. for different 
values of the applied current /. The quadratic surface i = is shown, and 
contains the cubic curve Sq = S a U S r , the family of equilibria of the associated 
fast subsystem, which undergoes a Hopf bifurcation at the point marked y = y sp . 
A canard cycle without head is shown, and follows So- The cycle grows into 
a maximal canard, developing a 'lip' as it passes over the fold of the quadratic 
surface. When the cycle grows beyond the Hopf point it leaves the surface and 
grows a spike that winds around So via the fast subsystem. 

We now demonstrate how the analysis of previous sections can be used, by 
applying the exponential microscope followed by pinching, to derive the birth of the 
first spike. For a general analysis of the spike-adding phenomenon in terms of slow 
manifolds of saddle type and the relative position of their (un) stable manifolds, see 

Non-spiking cycles remain close to the parabolic surface z = 0, so we begin by 
considering only the dynamics therein. (A stronger case for this is when the fast 
subsystem itself has two time-scales, so that (|60l) read e\x = y = E2Z = 
with ei < £2 < 1, and we can restrict to the fast subsystem by setting e-i — 0. This 
subsystem itself has slow-fast time-scales with ratio e = exje^ <C 1. In effect, the small 
parameter e pushes the Hopf bifurcation point in the fast subsystem in the direction of 
decreasing z. Recalling that the Hopf point initiates spiking, this causes the periodic 
attractor of the full system to lose its burst, following instead a small oscillation via 
the fast subsystem; such a periodic attractor resembles a canard with head where the 
orbit leaves the vicinity of the repelling branch of the cubic nullcline and relaxes onto 
the attracting branch by oscillating around it.) 

Let us then begin by fixing to the surface y — c — dx 2 , upon which the dynamics 
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is given by 

ex = c + I - y - ax 3 + (6- d) 2 , (61a) 

y =s(x-xi)-y. (61b) 

As for a planar system, we then straighten out the cubic curve where x = by 
replacing y with the coordinate h = ex, obtaining 

ex = h, (62a) 

h =7]d{x) + h( - Vd (x)-l) , (62b) 



where 

r) X (x) = c + I - ax 3 + (b - \)x 2 . (63) 

The nullcline h = bifurcates when = h = dh/dx — dh/dh, with solution 
= n d (x) = rj' d (x) — e = es + 2h(b — d) — 2he/x, giving 

x h = ^-(l±b-d± - d) 2 - 3ae) , (64a) 

h h = T£s/2\/(b-d) 2 -3ae, (64b) 
I h = ax\ - (b - d)x\ - I. (64c) 

At the parameter I — I hj a Hopf bifurcation takes place (as is confirmed by a simple 
stability analysis). The maximal canard closely follows the curve given by the nullcline 
h = when I = I h , that is h = ej (x), where 

10(X) - e-iJxY (65) 

In the previous two sections we then magnified the strip \h\ < e. Let us omit this 
(as we found it does not distinguish different canards from each other), and instead 
magnify around the maximal canard by introducing the variable W = (h — ejo (a;))' 6 ' , 
giving the system 

x = -W [1/£] +7o(a0, (66a) 

W = m 1 -^- ([ V ' d (x) - e - e^(x)][W^ + e^ (x)} + er, d (x)) . (66b) 

The nullcline bifurcates when = W = dW/dx = dW/dW, the solution of which 
gives 

I Q = ax 3 - (b - d)x% - I, (67) 
where x is a solution of 

n' d {x ) - £7o0o) = £, (68) 
and, for completeness, we have also 



Wo = „. dV s - £ 7o(x ) • (69) 

Beyond this value of I , the canard develops a head that remains, for a short range 
of parameters, on the parabolic surface z = 0. We can estimate this parameter range 
as follows. Let 51 = I — Ih, then the canards without head exist over the parameter 
range I — SI < I < Iq. The canard with head can exist over a similar range the other 
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side of Iq , that is, Iq < I < Iq + SI, during which the head grows over the region 

Vi < V < where y\ = I + c and yi = c + I + (|j (^p) ■ The head does not 
reach full size, however, being interrupted by the Hopf point of the fast subsystem and 
developing a spike. Let us say that this takes place at y — y sp , shown in hgure [T51 
If each spike grows linearly with I up to yi, then spike growth takes place over a 
parameter range I sp < I < I sp + ^ SI, and is preceded by a 'lip' — the canard 
briefly growing a head — over a range Iq < I < I sp . If we assume that the growth of 
the head with / is approximately linear, then 

I SP « Io + ^^51. (70) 
2/2 - 2/1 

It remains to find the value y = y sp at which the fast subsystem undergoes a 
Hopf bifurcation by solving x = y = and dx/dx + dy/dy = 0, giving y sp = 
I + c — axl p + (b — d)x 2 sp , where x sp = (b — yjb 2 — 3a) /3a. 

Finally, it is possible to capture the creation of the first spike via the canard 
explosion, by using an exponential microscope and looking at bifurcations of level sets 
associated with the flow. A complete description of the spike-adding phenomenon 
in such bursters, the necessary exponential rescalings, and the bifurcations of level 
sets, is beyond the primarily expository scope of this paper; these will be addressed 
in future work. We turn now to a more general result, a local classification of the 
different forms of canards that are possible in systems of arbitrary dimension. 



5. Classification of canards in piecewise-smooth systems 

Recently, a class of discontinuity-induced bifurcations affecting periodic orbits have 
been classified in piecewise-smooth systems. In [18], it was stated that these 
"catastrophic sliding bifurcations" included a classification of three generically 
occuring canards. We will show here how these can be derived from singularly 
perturbed systems by the method of pinching. 

In a neighbourhood of one of the folds {f'(x) = 0} of the critical manifold Sq 
in we can neglect terms cubic in (x — 1) and, after a simple change of variables, 
write 

1 2 , , 

£x = V ~ 2 X ' ^ la > 
y =q-x. (71b) 

This expresses a slow- fast system with a critical manifold So given by y = \x 2 , in the 
neighbourhood of a fold point at the origin. We can use pinching to derive a piecewise- 
smooth approximation to this system. Let h — y — \x 2 , so that the critical manifold is 
h = 0, and choose a neighbourhood, \h\ < a for some a > 0, to form a pinch zone. The 
nullcline h = 0, labeled S in section[3j is given by S = {(x, ^x 2 + e — lj) : a;GlR}. 

We can pinch in the y-direction along fibres (fTT)) given by Pa;(A) = (x, \x 2 + aX). 
This is equivalent to defining a new coordinate y = y — trsgn(/i) which, by substitution 
into (THI) . satisfies the piecewise-smooth system 

x=°sgn(y-^x 2 ) + 0(y,x 2 ) , (72a) 

V = q-x. (72b) 

From (TIT)|) . we can derive a sliding vector field F$ — 2-^(1,2;) on the switching 
manifold y = \x 2 . This appears to be divergent at x = 0, but from section [3l it 
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only applies on the sliding region, for which | ^ — 1 1 < £ , hence the point x = 
does not lie in this region so the divergence is never encountered. The boundaries 
are the tangencies of the upper vector field to the switching manifold at T 1; where 
x = q/(l + cr/e), and of the lower vector field at T 2 , where x = q/(l — cr/e), where 
the vector field curves respectively towards and away from the manifold. There is a 
sliding equilibrium (a point where Fs = 0) at x = q. 

If we plot the phase portrait in the {x, h) coordinates for varying q we arrive at 
figure fTBT a). Consider if the orbit marked ★ (for q < 0) is the local segment of a 
periodic orbit that detaches from the manifold at T%. Let q increase: the tangencies 
collide when q = 0, at which point the orbit can pass from the stable to unstable 
sliding regions; we call this a simple canard. For q > the orbit terminates at the 
sliding equilibrium, hence the periodic orbit has vanished, but has done so instantly 
(or catastrophically). 




simple canard visible canard 



Figure 16. Types of canard in two dimensions, as the parameter q varies. A 
small change of inset causes a jump in the outset. We classify these as: (a) simple 
canard, (b) visible canard. 

Notice that the tangencies T\ and T 2 in figure fTBTa) are such that the vector 
fields curve locally downward, confirmed by verifying that ft, < in the flows of (fTTj) . 
More generally in a piecewise-smooth system, the sign of h can be positive at both 
tangencies, or different at each, so the vector field can curve towards or away from the 
manifold. When all possible scenarios are considered, it is easily verified (provided 
h 0, see [TJJ or [TB]) that there are two other cases in which an orbit passes from 
stable sliding to unstable. One of these is uninteresting, because the flows either side 
of the manifold both curve towards it, keeping the flow within a neighbourhood of 
the manifold. The interesting case is shown in figure ITST b) . Here, an orbit with a 
tangency to the manifold is visible on both sides of the canard parameter, hence we 
call this a visible canard. 

The double tangency point through which the canard trajectory passes is not 
generic in two dimensions, and is only seen under variation of a parameter (such as q). 
In three dimensions, however, tangencies occur along curves, and generically these can 
cross to form double tangencies. We now show how these are derived from a singular 
perturbation problem. 

Consider the three dimensional system: 

ex = y-^x 2 , (73a) 

y = bz + cx, (73b) 

z —a, (73c) 

which is a local normal form for a three dimensional system with a slow variable 
added to (|7Tj) 25J. As above, we can let h = y— ^x 2 , define a pinch zone \h\ < er, and 



Canards and curvature: nonsmooth approximation by pinching 



26 



pinch along the y-direction by defining a coordinate y — y — asga(h), which satisfies 



a piecewise-smooth system 

x = ^sgn(y-^x 2 ) + 0(y,x 2 ) , (74a) 

y = bz + cx, (74b) 

z =a. (74c) 
The sliding vector field Fg is given on y = ^x 2 by 

x — (bz + cx)/x (75a) 

z = a. (75b) 



This can be recognised from [35] as similar to the slow subsystem of (|75|) when e = 0. 
Notice, however, that this piecewise-smooth approximation has been obtained for 
nonzero e. The sliding region is given by \{bz + cx)/x\ < cr/e, which excludes the 
line x — where (|75a|) is singular. The boundaries of the sliding region are where 
the vector fields in (|74[) are tangent to the switching manifold y = 0. They lie along 
z = (±j — c)x/b, labeled T\ and in figure [T71 These lines cross at the origin, and 
the sliding vector field is not well defined there. 




(folded saddle) (folded node) 



Figure 17. The 3 canards in three dimensions. A small change of inset causes a 
jump in the outset. We classify these as: (a) simple canard, (b) robust canard, (c) 
visible canard. Their corresponding classification in singularly perturbed systems 
1251 is indicated in brackets, noting that (c) seems not to have appeared previously 
in the singular perturbation literature. The vector field is tangent to the upper 
side of the switching manifold along T\ , and to the lower side along T2 ■ An extra 
arrow in (a)-(b) shows that the vector field curves towards the manifold around 
T 2 . 



We can infer the phase portrait of the sliding vector field by considering the 
desingularised system (xx,xz), which has at the origin: a saddle if ab > 0, a focus if 
ab < ab + jc 2 < 0, and a node if ab < < ab + \c 2 \ the focus/node is attracting if 
c < and repelling otherwise. We need then to observe that the equilibrium of the 
desingularized system (xx, xz) is not a true sliding equilibrium (that is F5 is not zero, 
and not even uniquely defined, at the origin). Note that the rr-factor reverses flow 
direction in the unstable sliding region in the desingularization, and that orbits in F5 
can reach the origin in finite time. These facts are usually indicated by adding the term 
"folded" to denote the origin as a "folded-node", "folded-saddle", etc. (degenerate 
cases arise at the marginal parameter values ab + jc 2 = and abc — 0, but these are 
nongeneric); see [25, 29]. Thus we arrive at two generalizations of the simple canard 
(figure fTBt a)) in three dimensions, shown in figure IT7T a)-(b). Both cases require that 



the Jacobian matrix 




of the desingularised sliding vector field (xx, xz) has 
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two eigenvectors in the sliding region, which is satisfied if c 2 + Aab > and 

a/e > (c+ hab/ (c± \J c 2 + 4a&)}) 2 , (76) 

which is found by requiring that the righthand sides of (|74|) . for the different signs 
in (|74al) . both point towards or both away from the switching manifold h = 0. The 
simple canard, figure UTT a). arises from the saddle-like sliding vector field (ab > 0). 
A robust case, figure fTTTb 1 ). arises from the attractive node- like sliding vector field 
(c < 0, ab < 0, c 2 + 4ab > 0), and is so-named because a whole one-parameter family 
of canard orbits exist. 

These were classified from piecewise-smooth systems theory in |18j . in which a 
third type of canard was introduced. Notice that, in figure HTTa)-(b). one vector field 
curves towards the switching manifold (along T2) while the other curves away from 
it (along 7\). This is the situation so long as (ab — a 2 /e 2 ) 2 > c 2 cr 2 /e 2 , and hence 
is always true if a/e is sufficiently small. More generally, however, the vector fields 
could both curve towards the manifold (if cafe ± (ab — a 2 /e 2 ) < 0) or both away 
from it (ccr/e ± (ab — a 2 je 2 ) > 0). Canards are possible in both cases, but only in 
the latter case can they depart the switching manifold in such a way that they can 
form one-parameter families of periodic orbits - these are the only ones of interest 
for canard explosions. The canard arises when the sliding vector field is saddle-like 
(ab > 0) with only one eigenvector in the sliding region f (|76[) is satisfied for only one 
of the signs ±). This we call the visible case (figure HTJc)), and generalises the visible 
canard from figure [TBT b). 

Figure [TTJ illustrates the phase portraits around the "sliding canards" that are 
generic in piecewise-smooth systems. Consider the orbits marked * to be local 
segments of periodic orbits, and allow their inset to change as a parameter varies 
(similar to figure [T6l but this time the local vector field does not change with the 
parameter). Then observe that when a periodic orbit deforms such that it passes 
through a double tangency point, it first forms a canard, and then, because its local 
curvature changes discontinuously, the periodic orbit will be destroyed. These belong 
to a larger class of four discontinuity-induced transitions called "catastrophic sliding 
bifurcations" [18] , the fourth of which we omit since it does not involve a canard. 

Figure [17] thus provides a classification of piecewise-smooth canards in three 
dimensions. Furthermore, in [Jj5] it is shown that these form a classification in K™ 
for n > 3. It is interesting to remark that the local sliding vector fields which give rise 
to this classification appear already in Filippov's seminal text on piecewise-smooth 
systems [11) . Figure [TrU provides their extension to planar systems, and in sections 
14.1114.31 we showed that piecewise-smooth systems, derived from slow-fast systems by 
the pinching method of section [3] typically exhibit the canard types classified above. 

A note of caution must be given on the unfortunate terminology that has 
arisen for this same singular point, in the different fields of singularly perturbed 
and discontinuous dynamical systems. The singularity lying at the origin of ([73]) . 
namely a non-hyperbolic point on a slow manifold, is called a fold (and in a addition 
a folded- node, folded-saddle, etc. depending on the slow dynamics), see e.g. [25]. The 
corresponding singularity lying at the origin of (|74p . namely a quadratic tangency of 
a piecewise smooth flow to both sides of a switching manifold, is known as a two-fold 
singularity, (this is because a quadratic tangency to only one side of the switching 
manifold is known as a fold), see e.g. [19]. 



Canards and curvature: nonsmooth approximation by pinching 



28 



6. Concluding remarks 

In this paper, we revisited the canard phenomenon in singularly perturbed dynamical 
systems, showing a close geometrical link between nonstandard analysis and piecewise- 
smooth methods. We first reviewed the exponential microscope approach of Bcnoit ct 
al. p], which allows the dynamics associated with a canard explosion to be resolved 
locally. We then showed that similar results can be obtained in the framework of 
piecewise-smooth dynamical systems using the concept of pinching. The exponential 
microscope renders the vector field almost vertical close to the critical manifold, 
providing an obvious choice of pinch zone, and we analyzed the parameter values 
where canards and maximal canards occur. We also applied these techniques to a 
three-dimensional burster, the Hindmarsh-Rose system, to analyse the transition from 
basic canard cycles to canards with one spike. The same mechanism can be found in 
other three-dimensional bursters (such as the Morris-Lecar equations [27] ) and has 
been recently studied from the standpoint of invariant manifolds, in particular, slow 
manifolds of saddle-type and their (un)stable manifolds [T3]. We have shown that 
pinching can be used to study the transition from canard cycles to spiking periodic 
attractors. 

The method of pinching introduced in section|31 and applied to canards in sections 
14.11 151 can be used more generally to approximate singularly perturbed systems by 
piecewise-smooth vector fields. In particular, we have used it to describe the dynamics 
around non-hyperbolic (fold) points of critical manifolds, where the Fenichel theory 
of normally hyperbolic invariant manifolds breaks down. This is an alternative to 
blow-up methods, which focus on regaining hyperbolicity by introducing an extended 
system. 

Pinching can be seen as the converse to regularization, where a piecewise-smooth 
system is smoothed out at the switching manifold. It has been shown only recently, in 
[26] , that the regularisation of a piecewise-smooth system is topologically equivalent to 
a singularly perturbed system, and that sliding orbits are homeomorphic to dynamics 
on a slow manifold. This result assumes normal hyperbolicity of the slow manifold, 
and therefore does not yet apply in the neighbourhood of non-hyperbolic points, such 
as those giving rise to the canards we are interested in (the resolution of this will again 
require blow-up methods). Ongoing and future work include the study of how well 
pinching can capture the full range of canard phenomena in R 3 , such as mixed- mode 
dynamics [24] [22], small amplitude oscillations near the fold [7], and spike-adding 
mechanisms [13] . 

Complementing this use of flow curvature to characterise canards, is the 
classification of canards derived from piecewise-smooth (or pinched) systems in 
figure 1171 A canard occurs where regions of stable sliding and unstable sliding meet 
at a point. The dynamics in a stable sliding region is non-unique in reverse time, 
while the dynamics in an unstable sliding region is non-unique in forward time (recall 
figure EJb)). When the two come together, orbits can be channeled from stable to 
unstable dynamical regions, creating the enormous sensitivity to initial conditions 
responsible for the suddenness of a canard explosion. 

Appendix: Proof of Lemma 5.1 

Given the system in ([34]) . 

ex = W [1/e] +£7o(x) 
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we differentiate along an orbit, 
dW 



^- = e(h-e 1 o) [£] - 1 (h , (x)-s 1 ' (x)) 
= eW 1 ^^ [h'{x) - £7o(a;)] 



where 



h'(x) = - = 



:(q -x) + h(l - x 2 ) 



= - h ( q -x) + l 



x h 

To find 7o(aj), note that £70 is a solution of (j34|) when q — q , so 

e-y' (x) = h'{x)\ q=qo = —(q - x) + 1 - x 2 . 

7o 

Substituting the difference between (fT9"|) and (|80j) into ([77)1. we have 



dx 



h 



(q-x) [q - a:) 



To 



and straightforward manipulation of the term in square brackets gives 



1 £ 
7o h 

1_ £ 

7o h 
- £7o 



(a; - 90) + 

q - 



7oh 



{x-qa)+e 

- qo 



To 



5o 



1 

To 



To 



+ £ 



h 

q-qp 



Finally, substituting this into ([8lj) we have the result of Lemma 5.1, 



W 



dW 
dx 



W 



To 



qo q-qo 

+ £ 



□ 



(77) 
(78) 

(79) 
(80) 

(81) 



(82) 



(83) 
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